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ABSTRACT 

Simulation techniques for nonlinearities in servomechanisms using 
the analog computer have long been subjects for investigations. In 
laboratory, these techniques play important roles for prediction of 
performance or in the study of nonlinearities in servomechanisms., 

Many successful methods in this field were developed already. 
Seven different basic nonlinearities have been chosen for investi- 
gation in this report. Net works for simulation of these nonlineari- 
ties are then generalized or initiated with complete analysis. Method 
of scaling is discussed, and illustrative simulation in a second order 
System was included to illustrate the ease of simulation with reasonable 
accuracy. 

The writer wishes to express his appreciation for the assistance 
and encouragement given him by Professor George Julius Thaler of the 


U. S. Naval Postgraduate School in this investigation. 


yop 





© YMBOLS 


© Output velocity 
S) Output position 
e Error signal 
% Diode forward resistance 
Ye Diode backward resistance 
£\ Amount of backlash 
(ws, Amount of coulomb friction 
D Amount of dead space 
t "pull in" voltage of relay 
a 1 iB 
2. drop out" voltage of relay 
S) Bar on top indicates computer value 
iV se 
Oi 
Mypt Maximum M for transient response 


M pw Maximum M for frequency response 
¥ System damping factor 
Op) 


G@) Subscript c denotes computer value 


Subscript p denotes problem vaiue 
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fr introadection. 

In most existing physical systems, it is usual to have some degree 
of nonlinearity in the system, however smail the degree may be. Or the 
system may de:iberately be designed nonlinear such as a relay system 
utilizing the available torque to it’s full advantage. 

If this is the case, equations describirg the existing physical 
system becomes nonlinear, the concept of roosts of a characteristic 
equation is not defined and there is no formal mathematical relation 


between the time domain and the frequency domain, since the principle 


of superposition does not apply, thus invalidating the basic tools for 


the analysis of linear system. For the transient response of the 
system, numbers to predict transient performance are not readily 
available. 

Thus a method of prediction of performance for nonlinear systems 
becomes very important in the field cf design of servomechanisms and 
to understand the effect of nonlinearities which may be present in 
the system; or in the design of nonlinear systems such as relay servo. 

Todate, the method of analysis of nonlinear systems is largely 
confined to phase plane method. But we are all aware of the laborious 
works and time involved with this technique as the order of differential 
equation describing the system becomes higher. And the error involved 
in the evaluation of transient performance from phase plane method in 
higher order system may become greater, 

One method of predicting system transient performance and fre- 
quency response directly at the laboratory is to simulate the problem 


to obtain solutions of equations that describe the system by use of 


analog computer. 





Although a number of basiz circuit to simulate nonlinearity have 
been introduced, a general tendency to avoid utilizing such circuits 
still remains. This report is dedicated tc encouraging the method of 
analog computer simulation in the study of the behavior of nonlinear 
mechanism which are commonly met in the usual servomechanism, by form- 
ing up the circuits, analyzing them in detail and by establishing scal- 
ing. 

These circuits are formed such that it will require the minimum 
possible number of amplifiers with reasonable accuracy for engineering 
purposes, Seven basic nonlinearities are selected for the subject of 
investigation in the report, end the method of simulation is based on 
the transfer function method to minimize the number of amplifiers and 
for ease of scaling for all nonlinearities selected. Section 2 con- 
tains separated circuits to simulate basic nonlinearities with de- 
tailed mathematical analysis and general scaling. Section 3 illus- 
trates actual usage of the cirucits discussed in section 2, ina 
second order system with resultant transient, phase plane and fre- 


quency response curves. The method of scaling is illustrated in 


sample simulation in section 3. 





Zs Circuits for simulation of basic nonlinearities, 

In this section, basic nonlinearities encountered in the usual 
servomechanism are represented by simulation circuits for use with 
analog computer. These circuits are analyzed and simulated in the 
Donner analog computer. Statistic curves were taken and attached to 
this report. 

Scaling in general....For scaling of a system with nonlinearity, 
it is best to neglect the presence of nonlinearity and calculate maxi- 
mum values of the desired quantities to establish a proper scaling for 
the system. In this way, the main circuit to represent the linear por- 
ticn of the system can be coupes after which the portion of the nonlinear 
circuit can be inserted. If the availability of components necessary for 
the nonlinear simulation circuit is limited, such as voltages of batteries 
or potentiometer size, it may become desirable to select scaling factors 
convenient for the formation of the nonlinear portion of the circuit. 

In that case, it would not require a great deal of time and effort to 
change scaling factors once established. 

The nonlinearity is then inserted in the linear circuit to form the 
main branch of the simulation circuit as in the case of backlash, re- 
lays and dead space etc., or added to the linear circuit to form a 
compensatory additional circuit of the system as in the case of cou- 
lomb friction. 

If the nonlinear circuit is inserted into the linear circuit to 
become a part of the main branch, the scaling of the nonlinear cir- 
cuit must be such that the scaling factors of the quantity be identical 


with that of the linear circuit at inlet and outlet stations. 





If the nonlinear circuit is added to the main linear circuit to 
become a compensatory additional circuit, scaling must be done in 
accordance with the equation describing the portion of the main branch 
to which the nonlinear circuit was added, 

These are covered in detail in section 3 of sample simulation 


of a second order nonlinear system. 


vo 





Ae Saturation. 
(a) Saturation phenomena in amplifiers etc., is represented by 
input vs. output relation as shown in (Fig. l-a). 
(b) Circuit to represent saturation. 


se Circule utiiizine batteriesmaisg. i-b). 


Where rdl, rd2 are diode resistances, which further are divided: 
rf....diode forward resistence, 


rb....diode backward resistance 


bad = + X>» =— Ao re — 
(Reo mye r -*Re EEO a, + CH TE ta. = ¢ 


me of =a a) 


—_— 


Ko = Kl, A = ra 





ue ae - fetes, } a a 
A ( Rt Ko Bn 7 i= Bo" lat Va, 
1 x, + E2 _ <r 
ein RK, Feo ~ fas 
ee = 3 = 
Ri 2 Yar 
case 1) 





case 2) XK. <- Ee. 








Vdey _ te 
-xX,=< _ kh” ee = ins 
rd, ee 


case 3) XM >t Ei 
rdiz=tbxoo rh=tero -x% 2 _ Ee, 





The result of this circuit is the negarive cof the desired relation 
eA 


shown in Fig. l-a. 


ii. Circuit utilizing power sources and potentiometers (Fig. l-c). 


Considering the upper diode while not conducting, the cathode 


potential "P" is; 


=, GE xcs yr ‘; a2 r ee) [Cina ant | 


me Caz.& 1 (-42)Xo J 
When P equais to zero, it conducts, 
Azx,E + (l—d2) Xp =o 
G2 J 
Xo = E 


I lp. 
Once the upper diode conducts, 


Az a) 
- -(%, ae eo 
I+ = I~ Az. 








Ro 
Provided the diode forward resistance is neglegible with respect 


to other resistive parameters. 


es R,-y alr , Re ar 


G27 a,v 
or Ry « | J a K | 
x, e- Se 


J~— Az 








And saturation voltages El, E2 will be: 
az 


Il—az 
(c) Static characteristics. 


a 





Static characteristic curves were taken with the circuit of 


= , = 1M. 
E, — ea — 4s Vv. 


Fig. 1l-b, where: R, 





Input vs. output plot. 

This plot was obtained by use of an K-Y plotter, 
with the input sine wave frequency of .01 cps. 
and the result is shown in Fig. 1l-d. The choice 
of input frequency was purely for conventence for 


the plotter. 








Fig. l-a Saturation characteristic 





Fig. 1-b Simulation circuit for saturation 





Fig. lec Simulation circuit for saturation 
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Be Dead space. 
(a) Free play or dead space about the reference point of the load 
is represented by input vs. output relations as shown in(Fig. 2-a.). 
(b) Circuit to represent dead space. 


i. © CircuLe nti lizinem batterieosm(ris. 2-b). 


let E.=E,= & Kes b 
\ 


case 1) JE, \Ee| >X, 


diodes do not conduct. 


. Xe 


a ® 


case 2) x <- cE. 
lower diode conducts. 


xX, - Re (-x + Ez) 
Jule | 


\ 
provided that the diode and battery resistances are negligible 


with respect to Rl and RO. 


case 3) , & | 
upper diode conducts. 
x ..& (X,- E.) 
0 Ri 
provided the diode and battery resistances are negligible compared 
to Rl and RO. 


The result is the negative of the desired relation as shown in 


Fig. Zeoen, 


ii. Circuit utilizing voltage sources and potentio- 


meters (Fig, 2-c). 
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Considering upper diode only, if the diode does not conduct, the 


cathode potential "P" is; 


xX 
(i~ar@ir a a 
P= = as oC (i-—a or - ar 


=~ G4E + (1-4)X, 
When P equals zero, the diode conducts. 
OE + Ci-adx, =o 
Gy 


= 6 a. 


im AA, 


while the diode is conducting. 











ae - ie 
xX, x ar C-adr b 
(i-ar tgp + & 
7 dr ar R, 
= 2 
if RK, i>~ a “> Se 
ats Ry 
Ay « ,* avy 
= Rk << | 
Xo = ~ Re XN \7 Ay 3 
and Ri C ae ) 
2.2. 
Ia, = 
Q 
1? ee 
*I£ slope b is defined as; 
bs Re , 
R, As ar 
a (=e a ee, 


Then the solution will be theoretically correct. 


/ XO Nae 
= oD R=? 


It was assumed 


that the '"a'''s were not corrected for potentiometer loading effects. 
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(c) Static characteristic curves, 


Static characteristic curves were taken with the circuit 


of Fig. 2-b, where; 


X= Jo 4mwt 


R=R, = | M. 


fi 


E,=EF,2en.5 WV. 


he inpub Vs, output plot, 
This plot was obtained by use of X-Y plotter. Input 
sine wave frequency of .01 cps. were chosen for pure convenience in 


the use of the plotter. Result is shown in Fig. 2-d. 
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Fig. 2-a Dead space characteristic 





Fig. 2-b Simulation circuit for dead space 





Fig, 2-c Simulation circuit for dead space 
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Seecauwlombp friction. 


(a) Coulomb friction at motor or load shaft is represented by 


the relation shown in Fig. 3-a. 


Where; Xl is velocity, 
XO is coulomb force. 


The coulomb friction can be referenced to the load or motor shaft 


as desired, and the simulation is basically similar, 
(b) Circuit to represent coulomb friction. 


iicmecult utilizine battermues (Pap. 3-b). 


Analysis of this circuit follows closely to that of the one shown 


mm Fig. I-c except the term ae goes to infinity in this case, and the 
\ 


outputs are shown below, 


case 1) xX, > O 
ee i 

case 2) pe tea 8, 
Xo = + & 


Where diode and battery resistances are neglected. The result is 


then, the negative of the desired relation of Fig. 3-a. 
li. Circuit utilizing power sources and potentiometers 


(ae. 3-C). 


Analysis of this circuit is similar to the circuit for saturation 


except the loop gain is infinity. 


case 1) xX, a 8, 


a Gr : 
Roe ian if —t <K | 


is. 











case 2) re > O 


AC ve 
af ee -¢ ee 
x, ee om E\ ¢ e K% 


j= Ay } 
Meow ofatic characteristics. 


Static characteristic curves were takea with circuit Fig. 3-b, 


X,= To an wt 
R,= 1.0 M. 
Sa eer 

i. Input vs. output plot. 


where : 


This plot was obtained by the use of X-y plotter, with 

input sine waveefrequency of .01 cps, and the result it 

shown on Fig. 3-d. 

(d) Discussion. 
This method of simulation requires velocity to be available. 

But this is the usual conditicn in simulation. Due to the open feed 
back loop at the instant when Xl equals to zero, the gain is infinite. 
The noise in amplifier will cause diodes to conduct one way or other, 
making output never to be zero. But this condition will not affect the 
actual simulation of coulomb friction because the physical system will 
behave in like manner and the coulomb torque applied in simulated cir- 


cuit while velocity is zero will not affect the output. 
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Fig. 3-a Characteristic of Coulomb friction 





Fig. 3-b Simulation circuit for Coulomb friction 





Fig, 3-c Simulation circuit for Coulomb friction 
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D. Simple backlash. 


(a) A simple backlash assumes two conditions, 


i. Output shaft does not drift when the gear is not in 


contact. 


ii. When the backlash is taken up and recombination occurs, 


output shaft velocity becomes instantaneously equal to 


input shaft velocity. 


These conditions automatically imply that the output shaft has no 
mass and has a great amount of friction in comparison with the inertia 
torque. 


The relation between input and output shaft position can be 
represented as in(Fig. 4-a.) 


Xl..Input position 
XO..Output position 


QA ..amount of backlash 


(b) Circuit to represent simple backlash. 


i. Circuit utilizing dry batteries (Fig. 4-b). 


case 1) YO+tE 
upper diode conducts 
Y= - (bx, + Yu) fe 
We y-t 
ee hy +E 


Provided the diode resistance is negligible. 


Ys _ 


~ Res 


4.|-2 





= - Rag Cb - Yu + E J “Lf -d 


19 





Y= —~bxi + E 41-4. 
l + ine € 
Substituting eq. 4.1-4 into eq. 4.1-2 
ae eee 
LL j-o 
=~ -bpRCSX, “+ RCS E 
ee ee 2 nS 
Ss Cob +E)Ci+Rc s) 
f + RCS 
¥.1-6 


bX, = Ee 


case 2) Y, <-E, 
Lower diode conducts. The analysis is similar to the case l, 


and it eventually yields; 


Xo = —bXx, + & 


case 3) Initially Yi = 0 


—--£<y¥ = 36x, Gee 


Yo =0 


diodes do not conduct, and 
oe Xo. O 


After diode has conducted; 


Y, #0 


20 


41-7 





and diode cuts off when; 
-E <~-Cbs +My) CTE 
~£ ie) | aa 
E = -(bx. t+ Yu) 
From eq. 4.1-4 
s(n iil 


l+Rcs” 
or 


Yu + ReosYy =-bX, + E 


“ by +Yys E- ROS Ya 4/—4 
Substituting eq. 4.1-9 into eq. 4.1-8 


“E=RCSYu-E 


ResM%=o 


Or, 2 (%)=0 


es Upper diode cuts off when Y4 reaches maximum. 
TE = —-PCEY, te 
RCS Yu =o 
A t 
aq 6 Yu)=O 


Lower diode cuts off when Y4 reaches maximum in the negative 


region. 


And if RC << | 


a x, = O 


ok = 
At when at Yao 


So, the above condition can be modified to say that when Xl 


reaches maximum the diodes cut off provided RC is very small in com- 


parison with unity. When diode cuts off Seo, 


ey 





Y; pow ~ Yy = 


= constant at it's maximum value due to charges 
in the feed back capacitor C, and X0 
value. 

slope b = gain constant of the first amplifier. 


backlash Am <FE 


ii. Circuit utilizing power sources and potentiometers, 


stead of the battery (Fig. 4-c). 


"i 7 
Upper diode conducts 
KY =~ (bX + YH) 


Ge 
f= A-a)p 


case 1) 


Ce 

atts a + \ 
(=a; 9 ar ' Re 
A % 00 


provided and diode resistance 


is negligible. 





nes Qe 4 Sr 
e+ a+ oe 
Let 
a a Oar ar 
= tit — — 
o ote R, + Re 
Then 
aE ~ 2e 
Ya" On- ee) 
—bxy 4+ = 
Nae i-a 


1b + JS RCS 
Substituting eq. 4.2-5 in eq. 4,2-4 


gf ae —_—— 
ie FON 7; E) 
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remains constant at the maximum 


in- 


M2 


Yn 4 


$¥2-3 


42-4 


4 .2-o 





Gey, > nears = ) 


~R¢S bx, + RCS E 
feet ics 


=-~(%+ Ya) 


ae 
- RcSbx + RicS aE ~bX + TA 
= ieee) ISS ag Ot ia 
l= Riess 
oe 
_ C14 RCS) be + =a EF) 
+ SRCS 
rE bx | (see discussion) , Xo mr -loX%, + = = 


case 2) Yi <-E 


Lower diode conducts, and the analysis is similar to case l, 


and it yields: 
A 
SS) COMO cepa | 
case 3) Initially Yy. = 


diodes do not conduct. 


en= 0 
Once a diode has conducted, 
Yy #0 
When a G 
— ae + ———- E 
=n ¢-6X1+ Yu) < i, 
diode cuts off. 


From eq. 4.2-5 
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h.2-6 





Y= bx + Fa E 
etd Kes 

K+ TRS = -bx + CE Ye-7 

I RCS Yy 


Q 
Yy + bX = a 


Comparing eq. 4.2-8 and eq. 4.2-6, 
or Y4 reaches maximum value, upper diode cuts off, 
if DRC KI 
dy so when *x%=0 
. + a3 


If Y4 reaches negative maximum, lower diode cuts off. 


4 Yu =o when X= 0 


gain censtant of #he first amplifier. 


backlash Qj = ao 


slope b = 


provided aad 


{c) Static characteristics. 


Static characteristic curves were taken with circuit shown 


in Fig. 4-c, where; 


YO ain Wt 
lo =| 
A= § Vv. 
© Biat 
Rx -0oFM 


24 





F= &s- 7 
y= 100K 
ae .is| 


_— 
= @ 2 @ = = 
— = 


(this value could not be made 
lower due to non availability of 
higher voltage battery.) 

Ll. Inpue@vs. output plot. 

This curve was obtained through the use of X-Y plotter. 
Input sine wave frequency of .0l cps. was chosen for the convenience 
of plotting, and the result is shown on Fig. 4-d. 
(d) Discussion. 

In any simulation of a servo system, the output velocity must 
be obtainable. In this backlash case, one has to have a differentiator 
just after the output Xo in order to obtain the output velocity. To 
avoid the erroneous result of the differentiator, one usually has to 
have some small resistance connected in series with the input capa- 
citor if the input is of the step nature. Once a resistance is connect- 
ed in the input side of the amplifier, the differentiated quantity will 
have error. If such is the case and the RCS in the denominator could 
be tolerated, the out put velocity is directly available in this 
simulation circuit. Y2 in Fig. 4-b and 4-c is the one, ie: 


From eq. 4.1-5 





Ys ~bRCeS X% A RcSE =O 
tr RCS 
if RO «| 


(while diode is conducting) 
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if — . 
ier = Qo (while diodes are cut of!) 


G fe rig fe aS ay a eS 
So, the volteze Y2 represents approximately the output veiccity. 


The quanity 3 


Q ar Cy 
&\ ar 
J 8 [este an Ry ar (<2 
CA OLY { Gy \ 
SS g es — «& 
Ot if j~ : [x 2 


(It has been assumed that the potentiometer loading effect is not 


corrected. ) 


The above condition can be realized with the use of higher vol- 


tage sources and lower resistance potenticmeters with respect to the 


value of Ri. 


Also, if a corrected potentiometer setting is chosen, and if 
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Fig. 4-a Characteristic of simple backlash 





r vii 


, Xo 


Fig. 4-b Simulation circuit for simple backlash 
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E. Ideal relay. 

(a) An ideal relay has no dead space and acts on an infinite- 
simal amount of the error signal. It can be represented by the re- 
lation shown in Fig. 5-9. 

Ales. <cberor ‘eitenal. 
XO....Relay output voltage 

This curve has exactly the same relation between input and out- 
put as coulomb friction. The only difference is in the quantity that 
input and output represent. For this reason the circuit used for simula- 
tion of the coulomb friction iS also uSed as a simulator of an ideal re- 
lay. 

(b) For circuit; analysis of it, and static curves, refer to para- 
graph C of coulomb friction. 

(c) Discussion. 

Due to the infinite gain in the open loop, the output voltage of 
this circuit will be El or E2 even if the input is zero, which seems 
to be contradictional to the characteristics of the ideal relay. 

In order to remedy this situation, reducing the loop gain has been 
tried. But then also, this will not be a true simulation of an ideal 
relay because of the finite slope in the output voltage to reach the 
desired level of El or E2 as illustrated in Fig. 5-6. 

If El and E2 are small, and if the gain is sufficiently large 
so that the finite slope is tolerable, there are three ways to obtain 
such a high gain. 


i. Decrease R1, (Refer to Fig. 5-c.) 
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This method seemed to be the best way to get high gain without in- 
troducing more noises among three methods. But Rl can not be decreased 
too low if voltage sources and potentiometers are used in place of 


battery as Shown in Fig. 5-d. 


hig — > | Say about 500; 


a OE an eS 


Ri Ia 
ar 
and R, may become not negligible. 


ii. Increase RO. 

Since the maximtm value of resistance available is 
limited, writer had to connect many resistors in series in the feed 
back path. This situation was rather an awkward one. It also tended 
to increase the noise effects. 

iii. Fractional output voltage feed back method. (Refer 
to Fig. 5-e). 

Xo ~ ~ Ro eX 

kK, a 
With the potentiometer coefficient '"a'’ very small, the gain can be 
boosted to a high value. But, for some reason, this method had more 
noise effect than either of the two methods discussed previously. 
Taking these factors in consideration, plus the fact that the 

physical system will also have noise such that the output will never 
become zero, writer has chosen the open loop circuit for simulation 


of an ideal relay. 


ou 





Xo 


Ey 
XK\ 


E2 


Fig. 5-a Charncteristic of ideal relay 





Vig, Seb Simulation circuit for ideal relay 





Fif. 5-<c Simulation circuit for ideal relay 
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Fig. Sed = Simulution circuit for ideal relay 





Fig. 5-e Simulation circuit for ideal relay 
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Relay with dead zone, 


F, 
Relay with dead zone but no hysteresis effect, has the input 


(a) 
output relation as shown in Fig. 6-4. 


VS. 
Xl....Relay input voltage. 


X0....Relay output voltage. 
D....Dead zone. 


(b) Circuit to represent the relation of Fig. 6-a. 


i. @ircuit using dry batteries. (Fig. 6-b), 
eme 1) —-&<¢ X, K+E 
K=O 
Xo= 
case 2) Da +E<KX 
Y= Xi7-E 


Koz -Ko¢x,-E) until ~ Re ¢xi- E) =p E. 
fy \ 


if XO becomes El; 
Ro Ke 


XxX, - Att (x) - E) —Ey 
\ 


if Ro te me 
Ko + V+ 
MX, = -& 
(> 
case 3) X%, <-2H =z 
if mak, + 


X> = Fe (x.-E) until fe (x, —E) = Ez 
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K; 
RB he 
Ro tk; ow 
= = f (X,\-E) + £4, ~ £2 
Ke 
hy 
rovided PPO, 
P h AS 
Result of this circuit will be the negative of the desired 
Gmawacteristics if Fo >? | 
| 
ii. Circuit using power sources and potentiometers. (Fig. 
6-c)., 
Az G3 [> 
X, <-E - — 
case 1) ae F 3 < < 3 Im Gs «te 
‘a | = . wip ae 
elie oe | a SE Gin, (i-ayk 
(i-as)rz u Q Y; rF — 
3% A 











emits Sale TAS 
IZ Fy 
Let 
igh Ri 
Then ‘ Gis 
=> -— Y= a3. ~ X- -— CE 
\ S23 ( \ map? I-az 3 
+ O03 % | 
X,% R(x 22 Bncit vo Pe Ea 
ow mm aa z)until oO ean. 
et IX. reaches = 7 
\~ ay. 
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ten XxX, = - Ko Gute ¢ yy, — 72 ED ae ss Ey 





ae eee, eos 
Be Oy =y, 
provided eae | 
a 
Sweet, Gale vey 
o ‘ 
case 2) E Az Az 
iE 
a I~ < nN < ° [az 


N=o Xe 


The result of this circuit will approximate the relations shown in 


Fig. 6-a in the negative sense. (2, does not affect the characteristics 


Smethis circuit if the condition Koy | is met.) 
R: 


(c) Static characteristics. 


Static characteristic curves were taken with the circuit of 


Fig. 6-b, where: 
K\ 70 am wt 
= =a 5S V- 
&, = £; S gas ve 


gain = Koo 


R, = | M 
Ros aod ff. 
feed back pot. setting a = .5 


i.) Input vs, output curve. 

This curve was obtained by using X-Y plotter with input 
Sine wave frequency of .01 cps. as shown on Fig. 6-d. 

From this plot, it can be seen that the diode cut off 
and conduction voltage has a little difference. If a person wishes to 
correct this effect, it can be done by uSing feed back as is shown in 
para. G of relay with hysseresis, 


(d) Discussion, 


Aside from errors arising from noncritical components of the 
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circuit, the most important factor in this circuit for successful 
operation with reasonable accuracy is to have the loop gain as great as 
can be obtained. But the boosting of this gain is limited, because of 
the dead zone associated with relay. If the gain is too high, simula- 
tion shows unstability in the dead zone due to noises amplified. If 
the gain is too low, it would not simulate the desired characteristics 
of this relay because of the finite time required for output to reach 
the relay output voltage. So, there is a compromise between these two 
factors. 

The degree of compromise vam depend on the computer and the 
problem itself. The gain was chosen as 400 in this case from the obser- 
vance of the performance of this circuit on the brush recorder. The 
output voltage jumped almost vertically to relay output voltage from 
zero without showing noises in the dead zone, with this gain on the 
brush recorder. The input voltage was chosen a low time rate of change 


voltage wave for this @bservance, it was .01 cps. sine wave of magnitude 


Ze Volts. 
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Fig. 6-a Characteristic of relay with dead zone 





Fig. 6-c Simulation circuit for relay with dead zone 
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G. Relay with dead zone and hysteresis. 
(a) This relay will have unequal "pull in" and "drop out" vol- 
tages, and it has the input vs. output relation shown on Fig. 7-a. 
X1....Relay input voltage. 
X0....Relay output voltage. 
- reullein” voltage . 


»+..' drop out" voltage. 


WIS No 


(b) Circuit to represent relay with dead zone and hysteresis. 
i. Circuit utilizing batteries (Fig. 7-b). 
case 1) +E >x, >-E 
Starting from zero, 
ye * Y= 0 
Xo» =O 
case 2) 
(1) When Xl increases from zero in the positive direction. 


o> +e 


Y, =—X 


Yee xr s-"+S 


where voltage drop in the diode is neglected. 


ie Be ¢x,-§) until me reaches El. 
o= R 


Then 


Xo # O 


, ~P 
Y =- Xi +b*o) xX, = & CX + bXo a, 


\ ies 
Xo = Recx,-£ ) 7.4-] 


Q \ 
1 7 8 








For the desired condition that the Xo jumps instantaneously to 
El or E2 from zero, gain Ee >>I as was discussed in previous para- 
( 
eraph, If the amplifier of the computer is a perfect one with no inter 
capacitance and a perfect 180 degree phase shift, the output Xo will be- 
have as indicated by eq. 7.1-1 during the fraction of a second when 
Xo is in the jumping stage to El or E2. 
But, due to the input and output diode limiter, current can not 
reverse in that same diode, and in the limit, eq. 7.1-1 will behave 
like eq. 7.1-2. oe 
Xp = ie (x- 4) | 
(until Xo reaches El or E2) 
But, an amplifier can not usually be perfect, then it is most 
probable that oscillation may occur during this stage, making it un- 
Stable. However unstable it may be, oscillation can not occur in 
this circuit because of two sets of limiting diodes, and this condi- 


tion will tend to make a steeper slope when Xo jumps to relay output 


voltage from zero. When Xo reaches El. 
K+ Ro 
ott 
Myon EEE (mw —£) + E, 


* 


\ 





where Y= diode forward resistance. 


«Kl, “tc 


\ 


Pp. 
Mo ~+FE at X= 2 =F while Xl is increasing. 


(2) After Xl has reached the maximum and decreases in the posi- 
tive region. ry Ry F 
XKo= Kot (x, +b&- 2) =e go E 
; 
2 6 = x a iL Ee 
4] 





Ye <e a input diode will cut off, 


X, tbe,<£=6 moc. . Se, 


we 


If this cut off voltage is to be A 





nies a g 
S-be = £ 
ee! 
bin oot 
Ey 
or TE, for x, > £ when Xl increases. 


ol 
; = when Xl decreases. 


case 3) x <-E 


Analysis is the same as case 2 except that the other diode con- 


duct and signs reversed. 


Xo=-E. for 


when increases, 
Xo = - guntil when decreases, 
ii. Circuit utilizing power sources and potentiometers. 


Griz. /-c). 


Starting from Xl = 0 and increases in the positive region. 








Xi=0 
ie 
for ~ 42 E,<xuc+ = F 
TZ=0 
Xp =0 
| 
G3 - 
when X E = 








Q A315 
(4+ SS te OL 
\- 3 ‘K 
eee |. pea 
Ye : ( mag J 
\ 
— Ro _ (x pee J 
X= El 6% - yD _ 


a Ey 
Bue. if S00, Xo will reach instantaneously to pom AY 
‘ 


and Xo will be limited at that value. The possibility of oscillation 


due to positive feed back is eliminated by diode limiters as discussed 


in the previous pages. 


When Xo reaches instantaneously to: 





Oe = > 
—_ = (EF 
As I= Ay 
Ro Qute L y C/3 = a oS ae 
xX, = Rot ast [ 5 + b ay * 7 a3 3J I Ay. 
K 
aE Oy ty é | Ay Vy << | v ~ Oy Ey. 
ll kOe , ay 


After Xl reaches the maximum and decreases in the positive region, 


a 
the diode in the input limiter will not cut off at K,= —=E, where it 
3 


started to conduct while Xl was increasing, due to the positive feed back; 


ie., a variable DC bias voltage. And it cuts off when; 














Ay Gz _ 
Se 2 acre = aan 3 
Or 
O3 OY. 
a Eo E 
may = oe 
if =e oe 
cuy. (? a | 
La =~ = 
2 7-2 @ 


This circuit will give similar results as the circuit appears 


on Fig. 7-b, provided following conditions are met. 
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(ce) Static characteristics. 

Static characteristic curve was taken with circuit appear- 
ing on Fig. 7-b with input sine wave of frequency .01l cps., where; 

a “ See 

EE, = &. =a. ¥ 

R= .1 M 

= 20 M. 

b = 3 different values. (setting of b was accomplished by 
setting potentiometer coefficient of Xo, while the gain 
of first amplifier was maintained at unity). 

Magnitude of input sine wave was chosen two different values to 
see the behavior of the circuit while Xo jumps from zero to El or E2, 
on the brush recorder. 

i. Input vs. output plot. 
This curve was plotted by X-Y plotter and is shown on 


Fig. /-e, where; 


f = ol aPS 


Fig. 7-3(1) 


Fig. 7-e(II) 
= eae f =: 


Fie. /-3¢11) 
fA... 


Zz” arm b= 


ie. a relay with dead space only. 
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(d) Discussion 
i. Since the nethod is based on using a varisble BC bias 
voltage, the success of the circuit will largely depend on the gain of 
tne second amplifier = 
Ki 

The greater the gain the better the simulation will be. But this 
gain is limited. Because, for higher gain, the noise amplification be- 
comes great and it appears in the output signal in the dead zone. 


inmeedas case gain 7 % 200 was chosen and the result was satis- 
} 


pa 


factory. However, there was a constant error in = ome indicated 
ome rig. /-e. 

To investigate this error, one more plot with b=0 (without feed 
back) was taken as appears on Fig. 7-e(III). As can be seen from that 
figure, it also has about the same amount of error, and the writer be- 
lieves that the error is most likely caused from the nonideality of 
diodes. 

ii. The error discussed above can be corrected. Fortunate- 
ly, the error is always outward from the origin, ie., the diode conduct- 
ing point is greater than the cut off point. Once this error is known 
to be present, it can be corrected by modifying the potentiometer co- 
efficient setting ''b". 

iii. A very low time rate of change voltage sine wave input 
was tested in this circuit to see the slope of Xo when it jumps to out- 
put voltage level, on the brush recorder. The result was a vertical slore 
even if the input was low time rate changing voltage. It is due to the 


unstability of this positive feed back network while the transition of Xo 


from zero to output voltage level. 
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Ro 
This steeper slope does not imply that the gain Ry is increased, 


and it will not give better results as a higher gain would. 

Supposing "'b'' is zero with the feed back loop left in place, the 
Giode must cut off exactly at the point where that diode has conducted. 
But, due to the finite gain = , there will be a finite time required 
for Xo to reach El or E2. In other words, there is a finite change re- 
quired in Xl in order to let Xo reach El or E2. The steeper slope of Xo 
does not mean that the gain is actually increased. [It only means that 
it shows a steeper slope inside the zone 3 as is described on Fig. 7-d, 


and it does not make the zone § narrower, which is the factor that 


governs the accuracy of the circuit performance. (Refer to Fig. 7-d). 
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Xo 


T= 
Xi 


3 


Fig, J-a Charactoristic of relay with dead zone and hysteresis 





Fig. 7-b Simulation circuit for rolay with dead zone and hysteresis 
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Vig, 7-c Simulation circuit for rolay with dead zone and hysteresis 





Fig. 7-d Period of unstability of the simulation circuit 
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am Sample simulation of nonlinearities in a simple second order 
S 


A simple second order linear system of camping factor -= is 
chosen to be simulated with each nonlinearity discussed in section 2. 

The method of scaling associated with each nonlinearity is illus- 
trated at each nonlinearity. 

Simulation of the linear system is also included in this report for 
comparison purposes and as a basic linear circuit to be simulated with 
nonlinearity. 

Method of simulation of the linear system used is the transfer 
function method although the direct component method also is adaptable. 

All scaling factors and circuit parameters in the transfer func- 
tion set up of the linear system are kept consistant throughout all 
Simulation of nonlinearities. 

In order for the above stated object, ie., to keep the transfer 
function of the linear system inside the closed loop to be consistant 
throughout all nonlinearities simulation, rearrangement in the basic 
system gain constant was made in the relay servo simulation. These are 
discussed in detail at paragraph F of section 3. 

Also, a,step input only was used in this report, and for each 
nonlinearity, three different plots are obtained. They are: 

i. Transient plot of output position and output velocity. 

In this plot, two different magnitude of step inputs are used with 
two different amounts of nonlinearity, and the recording was made 
through a brush recorder. 


ii. Phase plane plot of output position vs. output velocity. 
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In this plot, two different magnitude steps were used with two 
different sizes of nonlinearity, same as ones used for plot "a", 

X-Y plotter was used for plotting. 

i11. Frequency response curve. 

Two different amplitudes of sine waves were obtained from a func- 
tion generator and the signal was fed to the simulation board. For 
each signal, two different size nonlinearities were simulated. In 
this plot, input and output transient curves were obtained through 
the brush recorder from which the ratio of output to input magnitude was 
taken and the curve was plotted on semi log scale graph. The main 
object of the frequency response curve in this report is just to show 
that the Mpw changes in a system if the magnitude of input sine wave 
is varied due to the presence of nonlinearity in the system. Although 
a relation between the magnitude of input signal and the frequency 
response curve could be established in most cases, it is beyond the 


scope of this report and is omitted. 
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A. Linear system. 
) 


(a) Block diagrain (Fig. 1-1) 


(b) Mathematical analysis. 


6+76=KE 
Gc 26> ae We 
Bi St S17 225 7 6s*#+ ads + Ue 


W,, = 5 
ae 
% = 3 


Ue tp = v.42 6 
tp= aK- Xo 
M pt = ie pos 
For step input 0; = Y u(t) 


£6 mon) = ° Gal Sec. 
0, MA~AK — G64 T 


mee) Fee va 
Then 
tp = 2,22 Sac. Soe we 140" 


Tb won) = .927 (Winger = | 


(c) Scaling and circuit (Fig. 1-2). 
This circuit will be used for all other nonlinearities simu- 


lation without changing circuit parameters. 


Let: ra) 
i] 

Al 6, £ — 

me) 


a2 





mn 
any 
i 


ft 


ip 


And let: 


Oz = (.512)R34a=-532 Rzga00M Cra= |.04 mt 


Tile 
Gy = (-485) Ry Ct4 = 5085 Ry =1.00M. Coys i bat 


Mos =(.338)Re.3Cfz= 3468  Reog = 00 


(d) Discussion. 


For a linear system, the frequency response must be identical 


for all sizes of input sine waves. The frequency response curves 


shown on Fig. 1-5 verifies it with little departure. This deficiency 


was believed to have arisen from reading off the values on the brush 


recorded curves. 
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Frequency response data 
(linear system) 


ECcomPe) Wepre) Og M e. M 
OL 0943 21 E05 42 LS 
nO2 BSS 22. l,l 44 PelZ9 
-05 ~2oo 23,3 Th65 47 jel AP 5) 
04 aot 25345 1275 Sle3 baZe5 
nO ~472 28 1.4 ow aes, 1.44 
06 566 aye In 5 63 1.975 
07 . 66 s2.> 1.625 66 1.05 
O71 . 669 a2. 0 63 63.5 L655 
nO? 2 « ONS Cais, PeeL5 65 l .@25 
.08 ~754 225 1.495 60 eD 
ol 943 uo. 2 5 che 39 iS, 
5 1.885 3.4 iT 6518 aly 
4 S77, 8 04 i 0875 
1.0 9.43 Pe ial .0055 . 16 004 
6,=204m wt Oi = 40 am we 
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is 
SCSt1) 


K=2.2 5 


Fig. 1-1 Block diagram of the basic linear system to which 


various nonlinearities will be added, 





Fig, 1-2. Simulation circuit for the basic linear system, 
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B. Saturation. 

(a) ‘Block Diagram, 

Saturation in the amplifier can always be referenced to the error 
signal voltage, If this is done, the circuit to represent the satura- 
tion can be placed in the error channel as shown in Fig. 2-1. 

(b) Scaling. 

Let the amplifier saturate at: 


i, Weer 


ii, “= 


The output of the saturated quantity is still position error. So, 


the saturated quantity must be scaled the same as positional error. 





‘01 
(c) Simulation circuit (Fig. 2-2) 


Where: 


R. = ,> R= Ry = l.oo M. 


a Oa 
Gd, = 42= .335 
v. i 
E= ig 7G, MeO 
i Sai OB 3207 
jm yy 
vi 
\ 
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(d) Resultant curves, 

Fig. 2-3 and 2-4 shows transient responses of the system to 45 
volt and 60 volt step inputs respectively. The lower the Saturation 
voltage, the lower the overshoot and the lesser the output velocity. 
These are well shown in both transient and phase plane curves. 

If the input signal increases, saturation plays its part for a 
longer period and the above mentioned effects become pronounced, as 
is shown in Fig. 2-4. 

In the frequency reSponse curve, saturation decreases M pw by 
comparison with the linear system and the saturation effect is greater 
for a larger input signal in identical system, thus decreasing the Mpu). 
This situation conforms to the established theory. 

(e) Discussion. 

i. Error which may have been introduced due to the insertion of 
a saturation circuit will be due to the term a = -o8$33 ., 

Roughly this will give 3% of error in error signal. This error 
can be made very small either by using smaller resistance potentiometer 


or by chosing larger voltage sources. 
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(saturation) 
Véaty = roo: Ve 
6. M 
20 P00 
20.7 1.035 
22 ep 
24 Ws 
ads 1.35 
iD Wea e>) 
29°59 1.475 
30 Ue 2) 
30.5 as 3745. 
30 1b) 
16 8 
3.5 ae WB) 
5S) 025 
ell 005 


Frequency response data 
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Frequency response data 


(saturation) 
Vege) = 7O 
Q. M 
20 1.00 
20.7 1.035 
22.2 1.11 
24.5 1,225 
27 i305 
30 La 
31 £55 
sly 1.552 
Slee | P35 5, 
28.7 1.435 
18 ss 

2.0 13 
~45 0225 
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C. Dead space or Free play. 


(a) Block diagram. 


A dead space is assumed to be present about the zero position in the 
spring loaded load shaft, and the output position is measured at the 


load shaft (Fig. 3-1). 


(b) Scaling. 


The input and output of the dead space circuit are output position of 


the system, so that the scaling of dead space circuit must be done in 


accordance with the output position. 





DS D&D © 
KD AXE 
— D 
Dz A6 
Let 
(ia OS D = 3.32 V. 
;, Ome «| B= 6.667" 
(c) Simulation circuit (Fig. 3-2) 
Where: 
= 100 K 
R= woke 00M 
E=a45 v 5 a, 333 
G ; A= .035') = * 1a a 
ti gq aro6e4 a eemeci ae, SS Ce/ 
| a (= Gy ot: 


‘on 
i 
S 
ee 
= 
+h 





ee 

a, 2 oe 

6. 7 R, C A , ©, 
: 1 + &¢53 
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OC, 2 ~374 
(d) Resultant curves. -; 


Since the dead space was assumed to be present about the zero 
position ‘of spring loaded load shaft, it plays little effect on the 


system response for step inputs of the magnitude selected. 
The effect is almost indistinguishable on transient response curves, 


In the phase plane plot, it is better shown. If the size of dead 


zone becomes higher, output velocity at the start becomes higher. After 


some time has passed, the phase plane becomes almost same for different 


size of dead zone. 


However, it has a clear effect on the frequency response curves show- 
ing definitely different responses to different sizes of the input 


signal, thus indicating a definite presence of nonlinearity in the system. 


(e) Discussion, 


The slope of input vs. output of nonlinearity simulation circuit is: 


Re | 


Ri GAt Gar 
Ri 1+ Sf | 


as discussed on paragraph B of section 2, 
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However, the potentiometer loading was corrected in this case, then; 


slope = R x r 
ae Bias = ogee Or = ORE 
I-A e 
- ,9430 for a4 = ° 9684 
The % error is about three to seven %. But it can be reduced by adjust- 
ing Koso that Moy = slope. 
= ie RS! 3) ane 


Tt 





Frequency response data 
(dead space) 
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2 1.885 1.65 0825 4.5 »1125 
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1.0 9.43 04 002 - 006 0015 
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Frequency response data 


(dead space) 
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Fig e 30] 


Fig e 3-2 


Block diagram of the system with dead space. 
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Simulation circuit for the system shown in Fig. 3-1. 
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De Coulomb friction. 
(a) Block diagram (Fig. 4-1). 


The system open loop transfer function F(s) was chosen to be: 


a| 
Ris) = = os 


From this equation, the system equation without coulomb friction 





becomes: 
o. ee. ae 
If coulomb friction C is present at the motor shaft, equation 4,2 
must be rewritten to read: yb 
oe {- e Cc { 6 =K Ee ss 
Sis + > Oo nig oy Aaya — 
fa / 
Since il for the chosen system, and if 5 Q C _ , equation 4.3 can 


be rearranged to: 
/ : @ 
SH CstiI = KE ~C sAiqn Op 


where initial condition was assumed to be zero. 


Se, = KéE- Cane 4 Uf 
or = + | 
S6.=6, 


From equation 4.4, it is obvious that the transfer function method of 
simulation can be used for systems with coulomb friction if proper scal- 
ing is done. Then the block diagram can be established as shown in Fig. 
4-1, If the coulomb friction is present at the load shaft instead of 
motor shaft, coulomb friction C can be referenced to the motor shaft 

by letting cs Sic » where kg is gear ratio. 


(b) Scaling. 
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Let the referenced coulomb friction torque C' be; 


: C ee 
i | CD 
Where it is assumed that the unit is in accord. 
CAC = CC’ 
where a! 


/ 
Aa 
aC 6 = ome 
Cc 
7) 
One can choose any convenient value of C for more convenient values of 
; = 
AC - Usually C’is chosen for a convenient value of limiting voltage 
of the coulomb simulation circuit. Batteries may well best be used in 
the limiter, because of the scaling convenience. 


Going back to equation 4.4 in the problem form, 
6 = KeE-c4lgn ©, 
St! 


This equation can be written in terms of the computer values as equation 
ct ge a C ~~ 07 KC J 
fab -(——) ae [ke oe C yr 6, 
SP +t 

where, in this example: p( ¢ AAG 

From equation 4.5, the scaling procedure is clear and self explanatory. 


(c) Simulation circuit (Fig. 4-2) 


where: 


Let: od C’=.al2z 
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CG ' 12 jek, 
ll Cha (ge eee 
j= cfs 
(OF a= EtS — 6 O¢-/7 b 





Erio 
CE 
Ri shan. = Reve 5 
Rearranging a 4.5, 
= tl ewde - ~Txe atm OUTS J “1 
Let 4 
c! 
Ae = - = of Flt 
“" «bat 
Also, 
> We. h= (56 7 | R= 1.00 M. Ce = [04 wt 
2 


(d) Resultant curves. 

The effect of coulomb friction in the transient response of the 
system is as if the damping of system is increased. In Fig. 4-3 and 4-4, 
the overshoot becomes less when the coulomb friction is increased, while 
at the same time ouput velocity decreases. 

Due to coulomb friction in the system, the output position will have 


a finite steady state error when step input is imposed. This steady state 


error will be; OX Ex < & 

f 

: a/’= & Eos © UF 200s 
=o Ess § Yas V. 


This effect is well shown in the phase plane plot of Fig. 4-5. 

The frequency response curve shows a definite presence of nonlinear- 
ity in the system. 

(e) Discussion. 


In the case , aur. oolf » and the output of the limited quantity 


Ry. 
haste eLTOmor ot ce mes 
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This error could be reduced by the use of dry batteries in place of vol- 
tage sources and potentiometers, Also the same objective can be accomp- 
lished if lower resistance potentiometer or higher voltage sources were 


used, 
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Frequency response data 
(coulomb friction) 


C=) 
Ww 8, M 
0943 19.2 96 
1885 2022 1.01 
» 283 74) Sed 1.06 
o377 pe | 1.185 
472 26 1.3 
» 966 28 lee 
~613 26.5 1.425 
66 28 1.4 
, 724 25.8 1.29 
848 Zin 2 1.06 
© 943 16.5 0825 
1.885 3 Pa Ni. 
Be / oe 035 
5.66 26 013 
7.94 ~14 ~007 
9.43 005 0025 
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Frequency response data 
(coulomb friction) 


Cc’ = 107 
6, M 
18 9 
16.8 094 
20 1.0 
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24 le 
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26 Les 
25.8 £29 
Z3 1 
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3 5 
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025 20125 
14 .007 
05 0025 
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S, a 
5657 .918 
39 975 
42.2 1,055 
47.5 1.19 
52.5 1.342 
57 1.425 
D7 <2 1,44 
57 1.425 
3235 Jy J 4 
43 1.075 
34 84 
6.6 . 165 
1.5 OK Be) 
6 015 
0 24 006 
o15 .00375 





Gy x ; a2 Ss 


Fig. 4-1 Block diagram of the system with Coulomb friction, 
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Fig. 4-2 Simulation circuit for the system shown in Fig. 4-l, 
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Simple backlash. 


(a) Block diagram (Fig. 5-1) 


E. 
Backlash is assumed to be present in the load shaft of the system, 


and the output position is measured at the load shaft. 
The simple backlash implies that the load inertia is negligible 


compared to the drag force so that the load does not drift during 
load shaft velocity jumps tc 


backlash period and when gears recombine, 


that of the motor. 


(b) Scaling. 
The scaling of this problem is really simple. Since the quantity 


to be dealt is the position, and the circuit is to be inserted to the 


main branch, the backlash 4A must be scaled to the computer value by 


applying positional scaling factor. LA => 
o 
Let; bs 
ae = 3.33 
| “As -95 “> cole 
ae = =| 
=. = —— = 6.66 
(c) Simulation circuit (Fig. 5-2). 
Where ; 
EF =4¢5 
= ef M. 
C - | wf 
r P= GQ = 33° Gu oes] 
a 2 I~ <4 
I A ol cE. £667 A= 0688 
<= Ba 2 le 


From paragraph D of Section 2, 
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a 
_ CitecsyG 6+ Fa F) 


Oc aan & 4 Gi ar qi r 
1 t+ ORCS een, i. 3 


Since potentiometer loading effect due to Rl and R3 are taken care of; 
} = 
on -——- =, |». 0Y- for a= Oats | 


- of 
= |, oF for Ae. om 


Gye 6p Sue 


I-21 
kes sono OC, 2 034 wf Kee 5.02M 


6x bREC2 SS 
a adtd & = .395 
Ro OB 

(d) Resultant curves. 

it is a well established fact that the backlash has a destabiliz- 
ing effect on the system performance. From Fig. 4-3 and 4-4 of tran- 
sient response curves, it is seen that the overshoot and output vel- 
ocity becomes greater with more oscillations when the amount of back- 
lash is increased. 

Fig. 5-5 and 5-6 of phase plane plots, show more clearly of this 
situation, For the system damping factor of .333, of second order 
system, there should be no limit cycle, whereas the simulation shows a 
limit cycle for reasons that are not quite clear. 

However, with the understanding of the relation between maximum 
damping factor that will give limit cycle in the system, this circuit 
can be utilized for general engineering purposes. 


Frequency response plot on Fig. 5-7 and 5-8, shows that the in- 


crease of backlash corresponds to a greater Mpy . Mew is generally 
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higher than that of the linear system which is shown in Fig. 1-5. 

For larger input signal, backlash plays less effect on the system, 
thus becoming close to linear system performance, which in turn will in- 
dicate a decrease in Mow: These effects are well shown on the frequency 
response curves, 

(e) Discussion. 

i. If E in the Fig. 5-2 is increased to higher value, and if the 
potentiometer loading has been corrected, O can become closer to the 
unity. 

ii. As resultant response curves indicate, the sumulation shows 
limit cycle contradicting to the behavior of the physical system. 
Although the reasons are not quite clear, the writer suspects it to be 


caused from stored energy in the capacitor C of Fig. 5-2. 
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Frequency response data 
(simple backlash) 
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Frequency response data 
(simple backlash) 
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Fig. 5-1 


Block diagram of the system with simple backlash, 
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Fig.5-2 Simulation circuit for the system shown in Fig. 5-1, 
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F, Ideal relay. 

(a) Block diagram. 

The ideal relay is placed in the error channel as shown in Fig. 
6-1. 

In all previous cases, motor gain was assumed to be included in 
the constant K= 225, 

But, for simulations of a relay servo, it is assumed that K= 225 
does not include the motor gain constant, and it describes the open loop 
gain of the system not including motor. Also the motor is assumed to be 
an ideal torque motor. * 


The equations describing the system become: 


ee eo 2S 6: | 
ao, <Liso® E<O o% 


where T denotes the constant torque output of the motor. 


f 


= tae 


Bb 496 262-25 (7) E >O 
@+6 = 225-1 BO 


Then the transfer function set up used for simulation of all cases 
in previous paragraphs can be used without changing scaling factors and 
circuit parameters. 

(b) Scaling. 

Assuming that the torque T is in proper units the system equation 


can be arranged as follows; 


iu E70 
S 8, CitS) = 2-25 T 
6,= 2,25 T 


it S 
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ii. é<O0 


Ome 2 (6) 


i + S$ 


And the computer equations are; 


Ce 





6,6 = —2-25— 

= Ss + Pre At 
5 : on wo ase AT 
GAGs 225662 


aT fe Pe 
eT) £<° 


6- 3 


CG. 


From these equations, it is seen that the scaling is very simple, 


Similar to the coulomb friction case, 


If 


linear system simulation circuit can be used without rescaling. 


Let 


AT= XQ =-US , 


\ += 4S ae =: iS Cave 
| - 
7 P= 4S a = aren. 
(c) Simulation circuit (Fig. 6-2) 
where; 
kK, = 5 f2 M. 
laa Oo 
oe, OC 
\ tLe ad 
Gi. 
Gy Q,= -2 


(d) Resultant curves. 


Theoretically, the system must be zeroed in to the origin in the 
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phase plane. The simulated results shows limiting cycle contradictory to 
the theory. This is due to noise effect in the amplifier since the loop 
gain is infinity. 

In terms of transient response, the stability is poor. The stability 
becomes poorer if the motor torque is increased, exhibiting higher over- 
shoot and output velocity. These effects are shown in both transient 
response curves and phase plane plots. 

The frequency response curves are ploted as shown in Fig. 6-8 and 
6-9. The portion up to the resonance peak in the curve, is not a true 
response, because the output was an oscillatory wave. The nonlinearity 
is well verified to be present in the system by noting different res- 
ponses of the system to different sizes of input. 

(e) Discussion. 

Although the second order system with damping chosen must not 
exhibit limit cycle theoretically, the physical system will have noise 
involved, which in turn will drive the motor even when the error signal 
is zero. Considering this fact, the limit cycle observed in the simulat- 
ed circuit is not a detrimental factor. If a resistor is shunted across 
the feed back path, it eventually zeroes in to the origin, but not with- 
out a finite time sacrifice for the output to get up to relay output 


voltage. The writer preferred the infinite gain open loop circuit. 
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Frequency response data 
(ideal relay) 
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Frequency response data 
(ideal relay) 
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Fig, 6-1 Block diagram of tho rvstom with ideal relay. 
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Fig.6-2 Simulation circuit for the system shown in Fig. 6-1. 
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G. Relay with dead zone. 

(a) Block diagram (Fig. 7-1) 

The gain of amplifier Ka has only one effect on the system. It 
determines the sensitivity of relay dead space in terms of actual 
system position error. Thus, if the relay dead space is referred to 
a dead space in terms of error, Ka can be omitted. This implies that 
if the dead space simulation cirucit is placed before the amplifier 
with a proper scaling, the amplifier to simulate Ka can be omitted. 
Let ps width of dead space of relay in relay voltage. 

D & width of dead space of relay in terms of position. 


= 
= 


ap ey 


Wo=25 
The system equations become ; 
76g 16> T e>Ze 
J6+f6=-T ok 
JO+4O=>0 -—D<e< gz 
(b) Scaling. 
as Dy 
Ra 
Apteagz-ois 
== So vV- 
Since D is in terms of error = 
} o'=.-.4 D= 2045 D= oe 
Ye =r Ce Me ae OS ee 


(c) Simulation circuit (Fig. 7-2). 


Where ~— RP =\.00M 
yr = \ooKk 
E,= Fo Vv. 
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E-=200V 


Rae ro M 
QzHAivt 
gain of limiter amplifier = Re WES 
t 
bE2 => b,= -1304 








eo = 2 ee 1 a, . 08 613 for eS 
2 = [Ay - 
it o> © or 6 ae CO Te for O 
a = 4 
bit 
bir 
i - _ 
R, Ri 
lo, Ez 3 < E ee 
Xo a TT 


(d) Resultant curves. 


With dead space in the relay, system becomes stable. Wider che 
dead zone, the more stable system is obtained. These effects are 
shown in Fig. 7-3 and 7-4 of transientresponse curves, Because of 
noise amplified in amplifier, the limiter loop gain could not be gotten 
higher than 65 which is a little too small. Even with this gain, noises 


showed during the transition of dead space, Also in transient curves, 
drift of output position due to the noise effect is observed. 


In the phase plane plot of Fig. 7-5, two plots for different size 


of nonlinearity seems to fall on top of each other. This is due to 
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Similar isoclines of phase plane before and after the dead zone. As 
can be seen, two plots depart perceptively after a few cycles, 

In frequency response curves, the first part of the curve up to 
the resonance peak is not a true one because the output did not reach 
the steady state condition. 

(e) Discussion. 

Due to not high enough loop gain of 65, the simulation circuit 
imposes additional error inthe simulation. For error voltages less 
than e<Ze +2 region, the torque applied is not = 30 but 
(€- B) 65 . This error’can be reduced if the gain can be increased. 
If batteries or low noise voltage sources were used, gain may be gotten 
higher. The term BE sol imposes aebaer error in the circuit. If 

4 
© man = Go , the simulation will have 1.8% error in the torque. 
But this error factor can be reduced with the use of higher voltage 


sources or lower resistance potentiometer. 
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Frequency response data 
(relay with dead zone) 
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Frequency response data 
(relay with dead zone) 
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Fig. 7-1 Block diagram of a relay servo system with dead zone. 





Fig. 7-2 Simulation circuit for the system shown in Fig. 7-l,. 
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H, Relay with dead space and hysteresis. 
(a) Block diagram (Fig. 8-1). 
Let ; QO Relay pull in voltage. 


' & relay drop out voltage. 


k Z 

if 4 referred to error. 
a“ a 
/ 

Aa A referred to error. 
2 ~ 


The amplifier gain Ka has same effect as was discussed in the relay 


with dead space, and is treated in same manner. Then system equations 


become: 
E+ 6=0 e<ét lElincrease in + 
=e + {Oo a € r£ (Elincrease in + 
yorxté = O ae (El decrease in + 
JO +f6 =-T e<f (Glincrease in — 
pa:) +f © = @ E >¢g |\Eldecrease in — 
Let ; T=-67 
faz 20 
Li; , 
Feber Le ont T. 
z 
a % <a v Ase.ovs Vv. 
z= a 
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a a 
— 
Za 37. Z< ous 


(b) Scaling 


» 44 volt dry battery was available. 


Ge | 





fs 
=z zoo = ° 
a 7. 

A = ZAG ™ - 
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~ Pee. 

a -fqg=st 
ae 2 ye = | 

a ° 2ae& 

on Qh _ 

S-— =H. 


where; e 
Re -IM 
fig, Age 2ao™1 
E,=21°907 a 
d= -0o416 2 SEs £=s7. 


a a ~~ 
£-4-f f 
fem a -0454 
ap 
ii, eP_a 
. f=. OF for £-Loy 7 
K 
= me = 7 Oo" 
limiter loop gain = Re x 200 
\ 
(d) Resultant curves, 
eu 
The case of + = O; te, a relay with dead space only were 


included in this paragraph for comparison purposes, because i is 
different from the one used in paragraph F of this section, 


In the transient response curves,comparing Fig. 8-2 and 8-4 to 
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Fig. 8-5 (relay with dead zone only), it is evident, the relay with 
hysteresis has more oscillation with higher overshoot. This implies 
that the effect of hysteresis in relay servo is a destabilization. For 
cases of hysteresis chosen, it exhibits a definite limit cycle, while if 
£-deo, system is stable. 

Fig. 8-6 shows three cases of nonlinearities. If there is no hystere- 
sis effect in the relay, it zeros in to the origin as shown by curve 
“a'', while others reach limit cycle. 

The first part of frequency response curves up to the resonance peak 
is an approximation. The. output did not reach the steady state in that 
part of the frequency. The purpose of these curves are to show the 
presence of the nonlinearity in the system. 

(e) Discussion. 

As was stated in Section 2, the limiter gain z must be great in 
order for a satisfactory operation of this circuit. With rectified D-C 
power source in place of E2, gain could not be gotten higher than about 
50. It was then necessary to replace E2 with dry batteries to rid it 
of noises. Dry battery available in the laboratory was 44 volt. For 
this reason, it was found necessary to change the torque to T= 4 
volt from previous case of relay servo simulation. 

With the use of dry batteries, the gain could be increased to about 
200 which gave satisfactory results. 

In all phase plane plots, errors in "drop out" and "pull in" 
voltages is 12% in this simulation determined from the phase plane plots 
assuming that the plotter involves no error in itself. This error, is 


believed to be caused from inaccurate representation of 5 and the non- 
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ideal characteristics of diodes used. 

Also, the loop gain of non-infinity introduces error, In general, 
if the power sources are noise free so that the loop gain can be in- 
creased to a higher value, and if accurate representation of "pull in" 
voltages can be accomplished, it could be reduced. 

Or if the amount of error involved with the circuit is know}, then 


the error can be corrected by adjusting the potentiometer setting "f"'. 
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Frequency response data 


(relay with dead zone, hysteresis) 
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Frequency response data 
(relay with dead zone, hysteresis) 


T= 44 F=s 
£-2=47. 
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Fig. &-1 Block diagram of a relay servo system with hysteresis 
effects in the relay, 





voVys 
Fig. 2 Simulation circuit for the system shovm in Fig. 8&1, 
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4. Conclusion. 

Simulation technique of basic nonlinearity encountered in usual 
servomechanism is generalized and analyzed in detail, with illustra- 
tive sample simulation of second order nonlinear system. Utilization 
of this technique is apparently easier and time saving with reasonable 
accuracy for most engineering purposes. The scaling technique was 
illustrated, The error in general is less than 5% more added to the 
error involved in the simulation of linear system, except cases of relay 
with dead zone and hysteresis and simple backlash. But those exception- 
al cases can be made to have errors within the prementioned value. 

In all cases, error can be reduced by using more ideal components 
in simulation such as batteries, idealized diodes with more complicated 
setup cr non-ideal diodes, or potentiometers well balanced. It can 
also be minimized by more fortuitous choice of scaling factors as 
was discussed in each paragraph of section 3, such as making scaling 
factor of the time greater so that an error resulting from noncritical 
component of the circuit may become negligible in terms of the physical 
system values. 

The use of this technique in the study of the behavior of non- 
linearity in servomechanisms saves a great deal of time and labor 
especially if one deals with the transient performance and frequency 
response of a system, 

Given numbers for equations describing a physical system and 
numbers concerning nonlinearities, it is easier to set up circuit in the 
laboratory and simulate them to study or predict the system performance 


than drawing phase plane plots and tabulating data to further obtain 
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the transient response, When it comes to frequency response, the only 
mathematical approach "describing function method" does not give better 
accuracy, if as good. 

For reasons stated ete the writer believes that the simulation 
method in the study of nonlinear behavior is preferable to any other 
method, and this method concerning basic nonlinearities commonly met 


are studied and analyzed with illustrative scaling in this report. 
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